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In this thesis a numerical technique for solving 
the elasto-plastic torsion problem for Ramberg - Osgood 
type of materials has been prepared. The deformation 
theory of plasticity has been used to take into account 
the -work-hardening effect. The elastic and plastic tor-* 
sion equations have been satisfied in their respective 
regions and the location of the elasto-plastic boundary 
has been determined. 


Numerical results have been obtained for square, 
rectangular and thick L - shaped cross-sections. The 
effects of work-hardening on the elasto-plastic boundary 
and on the shear stress distribution over the cross - 
section, for different values of the angle of twist 
and the strain hardening parameter n, have been conside- 


red 



xii 


It has been shown that the elasto-plastic boun- 
dary is not appreciably affected by changes in the work 
hardening parameter n for low values of the angle of 
twist, but that the shear stress/distribution inside the 
plastic region keeps changing with the strain hardening 
pareraeter n. For large values of n (say 30 or more) so- 
lution approaches that for the elastic - ideally plastic 
case. 




1 . INTRODUCTION 


1.1 NATURE OF ELASTO-PLASIIC PROBLEMS 

In designing a structural member or a machine 
part, the engineer as a rule uses formulae which are 
based on the theory of elasticity. The highest stress ■ 
predicted by these formulae is then compared with the 
maximum stress which the material can stand without fai- 
lure. A factor, ; of safety is derived from this compari- 
sion. Now, in all but the very simplest cases the 
factor of safety obtained in this manner will differ 
considerably from the factor of safety defined as the 
ratio of the ultimate load to the design load. The rea- 
son for this is the fact that, the stresses set up by the 
ultimate load will exceed the elastic limit of the mate- 
rial 5 so that the theory of elasticity is unable to 
predict these stresses. In order to obtain the second 
factor of safety, the only one with a physical meaning, 
the mechanical behaviour of the material beyond the 
elastic limit must be ' taken into account . It is here 
only that the theory of plasticity plays an important 
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For most of materials it is not possible to 
describe the entire stress “ strain curve by means of a 
simple mathematical expression; so that in a given prob- 
lem the behaviour of the material has to be approximated 
by an idealized stress- strain curve, which emphasizes 
those aspects of the behaviour which are of importance in 
that particular problem. Some idealized models of mate- 
rial behaviour are described by rigid, elastic, perfectly 
plastic, rigid plastic, elastic perfectly plastic and 
elastic-plastic (strain hardening) materials'. These 
idealizations are useful for design work when the defor- 
mations are small and for carrying out detailed studies 
of the mechanism of fracture^ wear and friction of duc- 
tile materials under different types of loading. In a 
deformed body, when a part of the body has yielded, elas- 
tic-plastic models of stress- strain curves are used for 
calculating the stresses and strains. 

A complete solution of the general problem in 
plasticity involves a calculation of the stresses and the 
deformations in both, the elastic and plastic regions. 

In the former, the stress is directly connected with the 
total strain by means of Hooke's law. In the latter, 
there is no such correspondence and the stress-strain 
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differential relations have to be integrated by follow- 
ing the history of the deformation from the onset of 
plasticity at some pStint of the body. The process of 
plastic deformation has to be considered mathematically 
as a succession of small increments of strain, even when 
the over all distortion is so small that the change in 
external surfaces can be neglected. VThen the strains 
are large, the determination of the changing shape of 
the free .plastic surface necessitates the consideration 
of the deformation from moment to moment. , 

The state of strain in a plastically deformed 
element of a strain hardening material may depend on the 
entire history to which the element has been subjected. 
However, for the particular type of stress history known 
as proportional loading, in which the components of the 
stress tensor increase in constant ratio to each other, 
the strains obviously can be expressed in terms of the 
final state of stress*, for in this case the specification 
of the final state of stress also specifies the' stress 
history. The deformation theories of plasticity, which 
relate the total plastic strain to the final stress, are 
therefore applicable to proportional loading. 
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The solutions in the elastic and plastic re- 
gions are inter- related hy certain continuity conditions 
in the stress and displacements, which must he satisfied 
along, the elastic-plastic boundary* This boundary is 
itself one of the unknowns in the problem and is usually 
of such an awla^ard shape that even the stress distribu- 
tion in the elastic region can only be obtained by labo- 
rious numerical methods. The complete and accurate solu- 
tion of a plastic problem can only be obtained in a 
relatively few cases. 

1.2 REVIEW OF ELASTO- PLASTIC TORSION 

So far, all the work on elasto-plastic torsion 
has been limited to the elastic ideally plastic case 
only. The effect of work hardening ’ does not seem to have 
been considered by anyone. For an ideal plastic case 
(no work hardening) where the resultant shear stress at 
a point remains constant once it has yielded, an exten- 
sion of the membrane analogy was first suggested by 

-ri . . . . ■ 

A. Wadai (1923). This simulates the stress function 
for a partially plastic bar when the twist is so small 

that changes in the external contour can be neglected. 

^ ■ ' ■ 2 

iJsing this concept D.G. Christopher son (1940), and 

3 

D.G. Christopherson and R.V. Southwell (1938) found 
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the elasto-plastic boTindary for triangular and I-Sections 
by using relaxation methods. F.S. Shaw"^ (1944) calcu- 
lated the stress and the elasto-plastic boundary for 

5 

rectangular hollow sections. Sokolovsky;'- (1942) for- 
mulated an ingenious inverse method of solution for cer- 
tain oval contours. His method consisted in assuming a 
mathematically convenient shape for the elastic-plastic 
boundary and determining from this the approximate 
contour of the section. But as Prof. R, Hill has pointed 
out, the success of Sokolovsky'*'’^-' method was due to the 
fact that the elastic-plastic boundary happens to be an 
ellipse. Moreover, this method can be applied only when 
the plastic region has developed all over the boundary. 
Sokolovsky (1943), also investigated the plastic dis- 
tribution of stress in a conical bar and in a stepped 
bar consisting of two cylinders of different radii 5 but 
did not calculate the plastic boundary or the displace- 
ments. Eddy and Shaw (1949) applied relaxation methods 
to an approximate computation of the stress distribution 
and plastic boundary in a shaft with a collar. Apart 
from these, J.B, Walsh^ (1957) studied the elasto-plastic 

9 

torsion of a circumferentially notched bar and J.A.H.Hult 

(1957) derived the shape of the incipient region at the 

lO 

tip of the sharp notch in a twisted bar. K.R. Rushton" 
(1963) developed an electrical analogue solution of 
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elastic-plastic torsion for shafts having discontinuties. 
1.3 DESCRIPT lOI OF THE PROBLEM 

As pointed out earlier, all the previous work 
Was for elastic-ideally plastic materials only. However, 
most metals show a strain hardening effect when under- 
going plastic deformation. This happens because after 
minor shear slips have occured, they show no further plas- 
tic strains untill higher stresses are applied. Therefore, 
it is necessary to consider work hardening while analysing 
a plastic or elasto-plastic problem. 

To account for the work hardening effect, a 
deformation theory or flow theory of plasticity could be 
used. On theoretical grounds the flow theory is more 
sound. However, for small plastic deformations, the de- 
formation theory also gives good results. Further, it 
has been shown by B . Budianskey^^ (1959) that the defor- 
mation theories of plasticity may be used for loading 
paths other than proportional loading, without violating 
the general requirements for the physical soundness of 

the theory of plasticity. The' problem of plastic torsion 

12 - '■ ' ' ' ' ' ' 

was considered by J .H. HUTH (1955). However, his re- 
sults were inconclusive as they were based on a very 
coarse-mesh finite difference approximation. Later on, 
H.J. Greenberg^^ (i960) studied the comparision of flow 
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and deformation theories in the plastic torsion of a square 

14 

section for the Ramherg Osgood m-aterial. He established 
the validity of using the deformation theory in the tor- 
sion problem and the extent to which the predictions of 
the flow and deformation theories agree. Further, he jus- 
tified the use of the deform 3 ,tion theory by showing that 

the criterion due to Budiansky^^ is satisfied. He also 

14 

showed that for large va.lues of n, the Ramberg Osgood 
stress-strain law approaches the case of a perfect plas- 
tic material. 

In this thesis a method for solving the elasto- 
plastic torsion problem for Ramberg-0 sgood type of 
materials has been proposed, using the deformation theo- 
ry of plasticity to take account of the strain -hardening 
effect. For the plastic region of the cross-section a 
seperate torsion equation has been derived, which comes 
out to be a non-linear partial differential equation. 

Owing to the non linearity involved, it is not possible 
to solve this equation analytically. Therefore, a, nume- 
rical solution has been ■ attempted using the finite diffe- 
rence technique. It is an iterative procedure, wherein 
the elastic-plastic boundary and the numerical values of 
the stress function keep on changing after each iteration 
and finally converge to a stable boundary. In the final 
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solution the plastic and elastic torsion equations are 
satisfied in their respective regions. Further the values 
of the stress function match at the common boundary and 
also satisfy the condition of continuity of the resultant 
shear stress. 

The method has been tried for different types 
of cross-sections and for different values of the strain 
hardening parameter n. 

Since the basic nature of the plastic torsion 
equation remains the same for the piecewise linear work 
hardening model, the technique has been used for solving 
the problem of the, torsion of a square cross-section of 
a piecewise linear work-hardening material. The effects 
of the strain hardening parameter n on the elasto-plastic 
boundary and the resultant shearing stresses have been 
studied. The method has been tested for convergence to 
the elastic-ideally plastic solution by considering large 
values of n (say n =30). 

As the technique is based on the finite diffe- 
rence approximation, the effect of the mesh size has been 
studied from the point of view of computational time and 
convergence. This gives an approximate idea of the size 
of the mesh to be used for this method. 



All the numerical results, which have been obtained to 
illustrate the working of the technique, are based on 
the approximate properties- of steel. But the method 
could be used to study the behaviour of any other mate- 
rial in elasto-plastic torsion. 

All the results were obtained on an I.B.M. 
7044 Computer. 



2. basic theory of elastic ideally PLAgriC TORSIOH 


The plastic analysis of prismatic bars in tor- 
sion can most usefully be approached by first considering 
simple elastic torsion; as the mathematical formulation 
of the problem of plastic torsion (with work hardening) is 
similar to the elastic case except that a different stress- 
strain relation is used and a few more assumptions are 
made. Moreover, the plastic analysis of an elastic - 
perfectly plastic solid follows from the elastic analysis 
quite naturally and easily. Therefore, ST. VENAN’TS 
theory of elastic torsion has been briefly discussed below, 
before considering the general method of solution of elas- 
tic - ideally plastic torsion problems. 

2.1 ELASTIC TORSION 

A portion of a prismatic bar, shown in Fig. 

(l.a), is subjected to a twisting moment T. A system 
of rectangular coordinates x, y, z are used, where OZ 
is the axis of the bar parallel to the generators and 
Ox and Oy lie in the plane of the cross-section. Z = 0 
is one end of the section. Body forces due to the gravi- 
tational weight are neglected and it is further supposed 
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that local stress effects In the raaterial, iiamediately 
about the point of arplica.txpn of T, are redistributed 
in a short distance to give a stress distributJ.on across 
each section, uliich does not vary uitli Z. 

On deforra5..ng the bar each section rotates and 
also V7arps out of its initial plane. In Fig. (l.b), P 
is a point whose initial coordinates are (x, z) and 
which after to^isting moves f3?am P to P’ with coordinates 
(x + u, y + V, z + w) , u, V, w are small displacements. 
The angle of rotation of a cross - section at a distance 

z f rom t he origin is z©" , where e* is the angle of twist 
per unit length. For elastic defoliations, is small 
and is constant along the length of the bar. From Fig, 
(l.b) it follows that, 

X = r cos(S and y = r sin(2' ; so that 
Sx = - r sinp, dj^ , 6y ~ r cos |i dj?> , 
u ~ Sx ~ -y<?6 and v = 5y = xcSjj, 

Thus 

u - -y Z&- and v = x ze- (2,1.1) 

P’ is taken to have moved out of the original plane of 

P by an amount w. This amount varies from point to 
point, i.e. with x and y over the sect^ion, and can there,- 

fore be expressed as, 

w =®f (x, y) (2.1.3) 
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f (x, 7) is Called the wa.rpi.Dg function and is assumed 
to be independent of Zf so that warx)ing takes place at 
all sections} including the ends of the bar. 

From equations (2.1.1) and (2.1.2) the strain 
components can be expressed as 


£X = 

7>U- 

7>X 

= 0 

r 

= = 0 , €z = 

c>y 

0 

II 

Yxy = 

~dJ 

4* 

Y.V =: 

-2. X 

-Z©- + Z(9 =0 


Tyz = 

HZ 

4" 

"b W = 

■^y 

xe + 

^y 

(2, 1.3) 

Yxz = 

C) z 

4- 


- y©. + e- JX- 
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equilibrium for this ca.se reduce to 


+ "Mzx 


-^y 


X 


= 0 


( 2 . 1 . 6 ) 


Next, it is assumed that a function '^(x,y) 
exists which is capable of yi^ldiigthe shear stresses 
and satisfying the equilibrium equation (2,1.6). 

By tahing 

(2,1.7) 


^xz- and Tyz 




equation (2,1.6) is satisfied and equation (2,1.5) 
becomes 


-SxS ^y2 


- 206 


(2.1.8) 


or 


2 

V y 


- - 2Ge 


is ca.lled the stress function. 

By differentiating equation (2,1.4) it follows that 






l&S = G and 


2)y2 


On substituting these equations in eqn. (2,1.6), it follows 

2„ -^2. 


ITx +"dL, 

■Sx^ -5y2 


= v^f - 0 


(2.1.9) 


Now, supposing that the boundary of the section 
is defined by points (x,y) such that x = x(s) and y = yCe)^ 
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where s is arc length parameter measured from some fixed 
point on the boundary. In Fig. (l,c) a part of the sec- 
tion of the bar which is perpendicular to the axis of 
the bar has been shown. Considering the triangular ele- 
ment shown at the boundary, the shear stresses over an 
end face of the triangular area dA, are Tbcz and Czy. 
Together, these must be such that they produce no re- 
sultant force normal to ds, because there is no normal 
stress over the curved surface of the bar. This gives 


Tzx dx/ds 


( 2 . 1 . 10 ) 


Substituting from eqn. (2.1.7) in eqn. (2.1.10) it 
follows that 

^ =0 
"bx ds "by ds 

so that d^ _ pj 
ds "■ 

or V = constant, along the boundary of the cross 

section. (2.1.11) 

This is true for any boundary, whether it be an external 
one as in the case of a solid section or an internal one 
as for a hollow section. 
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The torque T, required to give a twist © will 
be given by, 


= 

X Tzy - 

y tTzx 

) dA 


=-((( 

+ 

y'dV 

. ) dx dy 


; ; 

ox 

J-y 





using eqn. 

( 2 . 

= 1 dy 

f ^ 

' dx 

~ f dx ( y 



J 15X 


J J 

'b7 


dy 

If the section is not a hollow one, then integration 
by parts gives, 

^2 

- Jdx -j^I|>dy 


T = - 


dy 


X ' 


74 


'Y “ ^ 

73 
(2.1,12) 

where, (x^, y^), (xg, y^), (xg, yg), (X 3 , y^) are the 
points on the boundary as shown in Fig. (l.d). 


Thus, 


Xf 


T = - 


dy 


Y dx 


- 

jsx - ygllfe - I Y dy] 


s^jfj Y dx dy 

lA. 


( 2 . 1 . 13 ) 


provided y , which is a constant along the boundary, is 



15 


taken to be zero on the boundary ( Vi ~ ^*2 ^ H^3 " ¥4 ” 

Thus for eqn, (2.1.13) to hold If must be taken to be 
zero on the boundary, 

Summai*izing the above results, if a function 

y (x, y) can be found which is zero on the boundary^!., of 

2 

the section being twisted, and satisfies "yif + 2 G© = 0 
then the shear stress distribution throughout the sec- 
tion and the torque T can be found. The magnitude of 
the resultant shear stress at any point will be 

2 


^ 2 

' y + 


( 




■& X 
grad y 


ClJL )■ 

?>y 


(2.1,14) 


2.2 ELASTIC - PERFECTLY plastic TORSION 

It was shown in the last section that |grady| 
represents the total shearing stress at a point. Phy- 
sically, this is the maximum slope at that point of the 
y surface (as in the membrane analogy). If K denotes 
the yield shear stress of the bar, then for yielding at 
a point 

jgrad'yj = K 
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Thus a limit is set on the value that the 
gradient at a point may attain. A further increase in 
the torque leads to a spread of the area over which yield- 
ing occurs. Constancy of the gradient in the membrane 
analogy may be simulated by erecting a constant slope 

roof over the section. An extension of the membrane 

15 

analogy was given by Nadai (1923 & *50). In this ana- 
logy it is assumed that the membrane under an increa- 
sing pressure is raised out of its plane, maintaining 
its largest gradient on the boundary and, having attain- 
ed the critical slope, proceeds to have increasing 
contact with the underside of the roof. 

In this manner a physical and intutive idea 
of the growth of the zones of plastic vielding under 
increasing torque may be had. From this analogy it be- 
comes clear that in general, yielding usually proceeds 
from the boundary inwards towards the centroidal axis. 

Throughout the outer plastic zones, the slope 
or yield stress is constant whereas within the section, 
the slope varies or the total shear stress decreases 
from K (assuming the material to be perfectly plastic) 
at the elastic - plastic boundary to zero on the axis 
of twist. To summarize, Ip satisfies, 



and 


-^2 ■tj; 7)2 ^ 

'o — “o = - 2G^ in the elastic region 

“dx*^ 'dj^ 

(^JiL)'+ ("^JSL) = in the plastic region 

'ey & 

and is zero on the boundary of the section. 

general method of solution 

For the convenience of mathematical discus- 
sion, the problem is better understood as consisting- of 
two parts. The first part consists of obtaining a 
solution to the elastic problem. 

Let R be a given region with p as its boundary 
as shown in Fig. (2. a). In R, it is required to solve 
the equation 

- “ 2Ge = - C (constant) (2.2.1) 

with 0 on the boundary p (2.2,2) 

Suppose that the answer to this problem has 
been obtained, then the quantity 

a = r (^ )^ + ( 1 ^ (2.2.3) 

^ J 

can be determined. From this it is possible to find 
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the maximum value of q in R + p . Then by multiply- 
ing Tj' throughout by ( K/q) the elastic solution corres- 
ponding to one point just reaching the yield point can 
be obtained. Thus, this part of the problem corresponds 
to obtaining the elastic solution over the entire cross- 
section 'When yielding has just started. 


The second part of the problem consists of 
considering the complete elastic problem (Fig. 2.b), 

The region R is now subdivided into two portions R^ and 
Rg (elastic and plastic), with external boundaries p 2 . 
and Pg . The position of the common boundary D, which 
is the elasto-plastic boundary is not known. Now in 
R^, equation 


2 2 




= - aC, a/ 1 


(2.2.5) 


has to be satisfied, where a is the ratio of the new 
angular twist q to <9c5 and C being the same constant as 
in eqn, (2.2,1) . ■ , 


On the boundary p' lii = 0 


( 2 . 2 . 6 ) 


In Rg, which is the plastic region, stresses 
will throughout be K only (material considered being 
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ideally plastic), 
this region is 

) 

On the hoandary 



Hence the eqn. to be satisfied in 
i 2^i 

4 - (^ )^ = K (2.2.7) 

j 

Pg , '^=0 ( 2 . 2 . 8 ) 


At all points on the cormon boimdary D, the 
requirement is ('^ )q = ( (2,2.9) 
But this condition is not sufficient to define the posi“ 
tion of D uniquely; a second condition is necessary. 

This condition is that at any point 0 on D 

( )q = (■^ )0 ( 2 . 2 . 10 ) 

“^1 c > 1 


Where, 1 is the common normal to the boundary. From 
eqns. (2.2.9) and (2.2.10) it follows that 






( 


~dVx 
"^y 0 


( 2 . 2 . 11 ) 


Now, eqn. (2.2.9) with either of the eqns. (2,2,11) is 
sufficient to define the position of the elasto-plastic 
boundary D. 
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The equation for the plastic region (eqn. 
(2.2.7)) merely states the fact that in (plastic 
region) the maximum slope of Ip is everywhere constant 
and has the value K. 

By virtue of equation (2.2,8) the contours of 
constant Ip will, therefore, be lines parallel to the 
boundary Thus equation (2,2.7) can easily be sa- 

tisfied for any region Rg and, infact, the solution is 
independent of the position of D. The portion Rg has 
been shown separately in Fig, (2.C) with its boundary p2* 
Contours of constant namely C^, Cg and Gg are shown 
at equal intervals d, which are all drawn parallel top2* 
On C^, Ip will have the value 

since ( 7 ) = 0 

- Kd (2.2.12) 

H being the measured normal to the boundary. 

The value of Ip on the other contours can be 
evaluated in a similar manner. In this way, working 
inwards from the boundary a solution of equation (2.2,7) 
for ip can be built up which will hold everywhere inside 
region Rg, irrespective of what happens in R^, and also 
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irrespective of the actual position of the bo-undary D. 
The solution depends on the shape of the boundary Pg 
and on the value of K, which is the yield stress in pure 
torsion. The real problem, ofcourse, is to find the 
elastic plastic boundary where the matching conditions 
given by eqns. (2,2.9) and (2.2,10) are to be satisfied. 

By adopting a iterative technique the second 
part of the problem can be solved completely, once the 
first part of the problem has already been solved. The 
entire discussion is based on a finite difference re- 
presentation which can be summarized as follows: 

(i) using the method which has just been described, 
a solution ip is constructed at the different 
mesh points. The two curves, shown in Fig. (3. a), 
represent the traces of the function ip' (x,y), 

(x,y) on a plane perpendicular to the x,y 
plane. A is /supposed to be the point on the 
boundary where the maximum value of q, i,e K 
occurs; and N is the normal to the boundary. 

The very method of constructing Ip makes the 
two curves tangential at A. 

(ii) Multiplying the elastic solution Ip of the first 
part of the problem by a known number a = 
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the first approximation to Ij? ^ is obtained as 
= a\|/ . This solution will have values 
greater than those of 1|; at some points; and the 
line of intersection of the two surfaces and 
1|; g_ can be determined by interpolation. This 
line has been indicated by point (l) in Fig.O.b). 
However, the condition that the two surfaces be 
tangential at the boundary is not satisfied. 

(iii) The position of the common boundary D is adjus- 
ted by keeping it on the "ij; surface and at the 
same time by keeping all the values of fixed, 
except those on the common boundary, until the 
tangential matching condition (2.2,li) is approxi- 
mately satisfied. This can be done by using the 
finite difference formula for first derivatives. 
In this way the boundary D shifts to a new posi- 
tion, as is shown in Fig. (3,b) by point (2), 

Civ) Now, since the position of D and the value of 
Ipg^ on D have both been modified, eqn. (2.2.5) 
is no longer satisfied at points such as (3) 

(shown in Fig. (3.b)), Keeping the new boundary 
values and the boundaries fixed, the solution to 
eqn. (2,2.5) is obtained in Rg by relaxation or 
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anyother method. This, however, gives rise to a 
violation of the "tangent^' condition. But, for 
this new solution the lack of agreement in the 
tangent condition will not be as bad as that in 
(ii) if the initial guess is good. 

(v) In this Way the position of the common boundary 
is readjusted as explained in (iii) and the com- 
plete procedure, described above, is repeated 
until both equations (2.2,5) and (2.2,11) are 
satisfied simultaneously. 

In this way a desired solution is obtained, 
since the requirements of equations (2.2.7), (2,2.8) 
and (2.2.9) are automatically satisfied. After a few 
iterations of the above procedure a complete solution 
of the problem is obtained. 

The solution of elastic - ideally plastic 
torsion has been thoroughly investigated by many authors 
for different shapes of the cross-section. The refe- 
rences for this have already been given in the section 
1. The detailed procedure which has been described 
above is due to F.S. ShaW^ . He used the relaxation te- 
chnique based on finite differences. 



foekulatioh of the torsion problem for 

WORK HARDENING MATERIALS 

The elastic - ideally plastic is a highly 

idealized model. In fact, there is always some work harden- 
ing associated with deformations beyond the elastic limit. 
Therefore, to obtain information about stress concentration 
and the plastic zone, particularly for the materials whose 
nature of deformation in the plastic range is far from being 
ideally plastic, strain hardening effects have to be taken 
into account. The primary purpose of this section is to 
consider the differential equation of plastic torsion. How- 
ever, since the coefficients in this equation are a function 
of the inelastic properties of the material, a brief dis- 
cussion of the ass'uifled yield criterion and non - linear 
stress - strain illations, as well as some comments on the 
deformation theory and strain hardening have been given, 
before deriving the differential equation governing plastic 
torsion. 

1* 

3.1 CRITERION FOR YIELDING 

A law defining the limit of elasticity under any 
possible combination of stresses is known as criterion for 
yielding or yield criterion. The two simplest and most wide- 
ly used yield criteria, which are consistent with experimental 
evidence, are the criterion of TRBSCA and the criterion due 
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to VOl MISES”', For most metals j Von Mises* law fits the 
experimental results more closely than Tresca's, but 
Tresca's criterion is simpler to use in theoretical appli- 
cations. In this thesis the Von Mises yield criterion has 
been utilized. For a general state of stress, in terms of 
the cartesian coordinates, this criterion postulates that 
yielding occurs, when 



( crx - (Ty)^ + ( (Ty -<rz)^ 


+ ( Vz - 


+ 6 




2 

Txz 


2 

't yz 



where, K is a parameter dependent on the amount of pre-strain. 

17 

As explained by Hill , certain assumptions and 
simplifications which are consistent with experimental evi- 
dence are made in arriving at a criterion of yielding. For 
the purpose of this thesis, the following assumptions are 
pertinent. 

(i) Time - dependent effects are not considered; 
this means, for example, that creep and 
thermal and viscous effects are dis-regarded. 

(ii) The Bauschinger effect is also disregarded 
since only a monotomically increasing torque 
has been considered. 
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(iii) The material is assumed to be isotropic, 
throughout its strain history, since large 
inelastic deformations are precluded. 

(iv) A hydrostatic state of stress has no effect 
on the yielding of the material. 

3.2 STRAIN HARDENING 

The concept of strain -hardening or work hardening 
in simple tension implies that the strain is a raonotomically 
increasing function of the stress. This means that the 
stress necessary for continued plastic deformation increases 
with the increase of plastic deformation. This coincides 
with the definition of a stable elastic - plastic material, 
Which requires positive work to be done during any cycle of 
loading. If beyond some point the material exhibits a fall- 
ing stress - strain curve then this portion would be classi- 
fied as unstable. 

A widely accepted general definition of strain - 

1'8 ■ 

hardening has been given by Drucker . Suppose a body is 
in equilibrium and subjected to surface fraction and body 
forces Now, by means of Sf^ and the stresses are 

first slowly increased and then slowly decreased by removing 
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an'd Sb^. Strain-hardening implies that for all such 
added sets of stresses the material will remain in equili- 
brium and that 

(a) Positive work is done by and Sb^ during 

the application of the added set of stresses* 

(b) The net work done due to 6 f.> 6b. during 

X ^ 

a cycle of application and removal is posi- 
tive if plastic deformation has occured in 
the cycle, or zero if only elastic changes of 
strains are produced. This definition of 
strain hardening is widely accepted since it 
is in agreement with experimental studies. 

Assuming that a strain hardening material obeys 
the Von Mises yield criterion and is the initial yield 
stress, then the radius of the Mises circle, or yield locus, 

J2 

is -g Yq. Further plastic itraining alters the current 
yield, locus. If the straining is continued to Y^ and the 
strains are then completely removed and if isotropy is assu- 
med throughout the cycle, the material now possesses a yield 

gfp 

locus which is a circle of radius -g Y^ . The latter circle 
will be concentric with the former circle. Such a material 
is called an isotropic strain - hardening material. The 
yield locus retains its shape but expands with the stress and 
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strain history. This concept of isotropic strain - harden- 
ing facilitates the mathematics of the problem, bait may in 
reality be considered only as a first approximation. 

3,3 DEFORMifflON THEORY 

The defoimation theory of plasticity is due to 
Hencky and Nadai. The term deformation theory was given by 
Il’iwshin. This theory assumes that, as in the theory of 
elasticity, there exists a one to one correspondence between 
the stress and strain. 

The popularity of the deformation theory may be 

attributed to its mathematical simplicity as well as due to 

its agreement with early experiments in which the stress 

ratios remained constant. Flow and deformation theories 

themselves agree if the stress ratios remain constant. It 

is well known that the defomation theory is inadequate in 

cases wherein there are large departures' from proportionality. 

However, when this does not occur, this theory can be con- 

18 

sidered to be less objectionable. Based on Drucker^s 
thermodynamic definition of work hardening it has been pro- 
ved by Budianskey^^ (1959), that the deformation theory is 
consistent with the flow theory for some restricted paths of 
loading apait frcm proportional loading. 
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A general stress strain relation, based on the 
deformation t]^eory, which was proposed by W. Prager^^ (l945) 
is given by 


where , 


£= f (Jg, Jg) 


p (Jgi ^3 T + q (Jg, 

L. 


4^ i 

(3.3.1) 



- 2/3 I 


(3.3.2) 


and p, q are polynomials, which are homogeneous in the stress 
components. 


Prager> has also discussed some special cases of 
the general relation given by eqn. (3.3,1). One of the 
simplest relations is obtained by putting p - 0, q = 1 and 
by assuming f to be a function of Jg only. The resulting 
relation is of the form, 

£= f (Jg) £ (3.3.3) 

This particulap stress strain relation stipulates that the 
ratios of the principal strains equals the ratios of the 
principal components of the stress deviate r. 

3 . 4 non-linear elast ic plast ic behaviour 

The non-linear elastic plastic behaviour, which has 
been considered here, is of the type discussed by Ramberg and 
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Osgood (1943) for explaining the tensile stress-strain 
relations for certain materials. Ramberg and Osgood have 
shovm, that most of the stress strain curves in tension can 
be represented by the relation 


¥■ 


1 + ( 


o 



(3.4.1) 


If all the stress components are non-dimen sionalized by the 
yield stress and the strainsare normalized with respect 
to, the corresponding yield strain, ~ > then eqn, 

(3.4.1) reduces to 

£= 0^ (3.4.2) 


where^ 0.02 (approx.), corresponding to the usual engi- 
neering definition of yield. 

Tensile curves for 0.02 have been shown in 

Fig. (7), for several values of the strain hardening co- 
efficient n. 

The application of eqn. (3.4.2) is restricted to 
monotomically increasing stress or more commonly to condi- 
tions of no unloading. 

In this thesis, only the simplest deformation theo 
ry has been employed. Pla'stic deformation is assumed to be 
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It 

independent of the hydrostatic component of stress 3 » 

and is also assumed to he completely determined by the second 
invariant of the stress de viator 



(3.4.3) 


(The first invariant being Identically zero) 

The invariant has been introduced in the form of the "effec- 
tive stress" defined by 




S.. 

13 


(3.4.4) 


Thus in simple tension <s^ - cr" 5 and the Mise’s yield crite- 
rion, for this theory, reduces to CJ^ = 1. Ofcourse, a 
smooth representation such as given by eqn. (3,4.2) admits 
only an effective or approximate yield condition. 


Now, a generalized stress-strain relation which 
reduces to eqn. (3,4.2) in simple tension is 




(3.4.5) 


This is the final form of the stress- strain relation which 

has been made use of in deriving the equations of plastic 

torsion. This type of a relation has recently been used by 
20 

many authors. 



PIECEWISE LIKEim PvELATION 


A somewhat less realistic representation of the 
tensile stress - strain "behaviour of common metals, other 
than the Ramberg - Osgood relation, is the piecewise linear 
relation. However, this approximation does model certain 
features of plastic flow. In particular, the behaviour of 
Alluminium alloys is very closely approximated by this type 
of an idealization. The assumption of piecewise linearity 
states that in each plastic regime (i.e side or corner), 
there exists a relationship, between the stress rates and 
the strain rs^es, whose coefficients are material constants 
which are independent of the, strain history of the mater- 
ial. In addition, the effect of strain hardening upon the 
behaviour of the yield condition itself is also linear* 

For such an idealization, the general stress - 
strain relation is given by 

^ 55. ^13 h 3 

where f - l / ^ 
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Using this relation, the ecjuations of plastic 
torsion have been derived for the piecewise linear beha- 
viour of materials also . 

3.5 plastic - TORSION 

In deriving the equation for the plastic tors3.on 
it has been assumed that the basic assumptions of the theo- 
ry of elastic torsion are still valid for the elastic - 
plastic state, except that a different stress - strain re- 
lation has been employed. That is to say that, even on 
the portions of the cross section which have become plastic 
the only non-zero components of the stress and strain are 

Txz > Tyz Ykz > Ttz • 

Further, in order to keep the mathematics of the 
stress - strain relations as simple as possible without 
losing the essential features of the theory, the plastic 
meterial has been considered to be incompressible. All 
viscosity effects have been neglected, that is, the stress- 
strain relation has been assumed to be independent of the 
speed of deformation. For moderate strain rates and tem- 
peratures both the assumptions are well 
structural' metals. 



Now, starting with equation (2,1.3), it is 


clear that 


■€ 


X 


^ ~ -£ = ^ = 0 
y z xz 


(3.5,1) 


and 




xz 


^w 

'bx 


ey 


(3.5.2) 


Xiz = 


Because only two components of the stress, namely 

T'xz Tyz "been assumed to he non - zero, the 

first invariant of stress, 

I = (TT; = 0 (3.5.3) 

kk 

Substituting zero for in eqn. (3.4,3) and (3,4.5) 

gives 




'^3 


(3.S.4) 


and 


"£13 




(3,5.5) 


Replacing by from eqn. (3.5.4), and expanding 

equation (3.4.4), the value of the effective stress comes 
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Once again, assuming tlxat a stress function exists 
such that 



1X2 




(T^ 


(3.5.7) 


the expression for ^ reduces to 

^ e 

The generalized stress - strain relfeion given by eqn. 
(3.4,5) gives 


- (1 + V ) T- + 1 '^ - 7 - 

T xz I xz (- xz 

and ~ q n-1 _ (3.5,9) 

£y,= (l.V)ry,- H 0^ Ty. 

Which on substituting the values of the stress components 
z ®in. (3.5*7), reduces to 


(1 + V ) + |c< 


% / (ff 






Next, by expressing the strain components in terms of the 
displacements and its derivatives from eqn. (3.5.2), eqns. 
(3,5.10) gives 
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-6y)E 


1^(1 +i) ) + I 


and 


y 

(3.5.11) 


n-1' 


1 ( + 0x) E = ~((1+D) + ^oC 

2 7)7 ] 2 " e j 5fx 

^ ^(3.5.12) 


<r^ V'V, 


Differentiating eqn. (3,5,ll) with respect to 
y and eqn. (3.5,12) with respect to x and subtracting, we 
get 

-0E = 


(1 +1) ) + ^ J y ^ (n - 1) 


5 ?“'" ( ( Yy ^ Y 


X 

J 


Or 


(3.5,13) 

- n-2 


~©E = l^d +1^ ) + f ^ y Y cr^ 

h ^iVx b ^ ^Vy V ■ ^ vV] 


(1 +V ) + I ^(2 - n) v- 7 ?. 


3 , n-2 

5<<,(n - 1) (T- 


^ <?~eVx^x 

(3.5.14) 


This is the final form of the differential equa- 
tion which is supposed to govern torsion in the plastic 
region. It has been put in this form, instead of writing 
in terms of stress function Ip and expanding all the 
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terms, to suit a special purpose which will he clear while 
discussing the numerical solution of this equation. 

Similarly, for a piecewise linear stress-strain 
relation a separate plastic torsion equation can he derived 
and it comes out to he, 


.0E 


(,rSVxh 


+ C <r^ )y.- 


where , 




and 


CTe 


3 
2 

= 0 
- i 


E’ 


- 1 


(3.5.16) 


if 


(7e > 




if erg ^ 1 
2' 


3.6 NON -DIMENSION ALIZAIION OF THE LENGTH SCALE 


The coordinates will he non-dime nsionalized as, 


- 2L- 


I, . \= f- 

where, L is a characterstic length of the prismatic har. 
Further a new stress function ^ defined hy 

^ ^ ,1^ ) = ^ ( X, y) (3.6,2) 

will he used. 
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Substituting for the variables in terms of the 
non-dimensional coordinates and , eqns. ( 2 . 1 . 8 ) and 
(3.5.4) for elastic and plastic torsion respectively 
reduce to. 


and 




'0E 


- 2G6 


(Elastic) 


(3.6,3) 


= Ijd + I (2 - n)^ 


2 

Vi. 




+ |<<(n ~ 1) 


The stress components are given by 


Txz 

and therefore 




yz 






(3.6.4) 


(3.6.5) 


is 


j ( 

(3.6,6) 

Equation (2.1,13), -which gives the value of the torque, 
becomes 


2 L’ 




d^ d)^ 


(3.6.7) 


Only these equations have been used in the solution of the 
problems and in the computation of the various quantities, 
such as the equivalent stress, torque etc. 



method of numbrigal solution 


There are many numerical methods for solving 
partial differential equations. Of all these, the method 
of finite-differences is used most widely because of the 
ease and universality of its application. This technique 
has proved its worth in practically every branch of science 
and technology; and the method has already been used by 
many authors for the solution of elasto-plastic problems, 
including that of elasto-plastic torsion. In this section 
of the thesis the method of finite differences has been 
briefly outlined and the finite difference equations for 
the differential equations (3,6.3) and (3.6.4) have been 
derived. Finally, the scheme for obtaining the elasto- 
plastic solution of our problem (torsion), has been dis- 
cussed in detail. 

4.1 METHOD OF FINITE - DIFFEEENCES 

The essential principle of the finite-differen- 
ce technique is to replace the differentials of a variable 
by the differences taken over finite intervals. The field 
of interest is, therefore, replaced by a net of grid lines. 
The points of intersection of the grid are termed the 
"nodes” and it is with these nodes as the focii that the 
differential equation is integrated. The finite difference 
technique, therefore, represents the 3?eplacement of the 



40 


differential equations by a system of difference equa~ 
tions which are essentially a set of simultaneous alge- 
braic equations involving the unknown variable at the 
nodes only. The solution of these simultaneous equations 
gives the numerical values of the unknown at these nodal 
points. The selection of the method for the solution of 
these equations depends on the order and the nature of the 
coefficient matrix of these equations. 

For writing the difference equation three types 
of difference operators are used; the backward difference, 
the forward difference and the central difference opera- 
tors. The central difference operator, being superior to 
the other two, has been used here; its order of convergence 
being h , (h being the grid spacing). The accuracy can be 
improved by using higher order finite differences but this 
requires the consideration of fictitious points beyond the 
boundary, for almost the same rate of convergence. 

The specific finite difference scheme for solving 
equations (3.6.3) and (3,6,4) has been described below. 

4.2 FINITE DIFFERENCE EQUifflONS 

ELASTIC TORSION EQUATION: Equation (3.6.3) which governs 
elastic torsion, is Poisson’s equation and has been solved 
by many authors, by using the finite difference technique. 
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The right hand side of eqn, (3.6.3) being a constant 2G© , 
for a given angle nf twist, we need only write in 

terms of finite differences* Salvadori^^ has given the 
finite difference operator for (Laplacian operator) 

in rectangular, polar, triangular and screw coordinates. 

¥e consider here only a rectangular grid system as shown 
below. I fl I 1 


I 


2 



5 -Rs 

O -fii 

■^4 

t 



4- 







Let, 0 be one of the typical nodes at which the 
differential equation has to be satisfied; and let, 1, 

2, 3 and 4 be the neighbouring nodes around 0. Then we 
know that (see Ref, 22). 
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rii +i3i 
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hi hg 
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hg+h4 


42.^4 

hg h4 


■S'" 


+ Error term 0 (h'^) 


(4,2 ,l) 


where, h^, h^, and h^ are the distances of the neigh- 
bouring nodes 1, 2, 3 and 4 respectively. 

Using this difference operator given by eqn. 

2 

(4,2,1) for "sy and neglecting the second order error term 
0 (h ), equation (3,6.3) can be replaced at the pivotal 
point 0, by the difference equation, 


-2 


\ ^3 




"^2^4 


2 


h^+hg 



•f -~:2 

V^4 




-* 2Ge 


(4.2.2) 


If the neighbouring points are at equal distances from the 
point 0, say h, the equation (4.2.2) further reduces to 

^ i’o ■ $1 ■ $2 “ $3 ■ $>4 " (4.2.3) 


Hence, the differential equation (3,6.3) has been 

transformed into a difference equation (4.2.3). Finite 

2 

difference operators for S7 can be chosen in other coor- 
dinate systems also and a difference equation similar to 
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(4.2.3) Can "be derived, 

plastic -TORSION EQUATIONS The core of our computational 
procedure is the derivation of the finite-difference 
equation for the plastic-torsion equation (3.6,4), To 
facilitate the derivation, the L.H.S. of the equation has 
been considered in two parts. The first part is 




Once again a rectangular mesh has been chosen 
(as shown on page 41) and the difference equation at the 
node point 0 has been obtained. Finite difference equa- 
tions have been derived separately for equations (4.2.4) 
and (4.2.5) which, when added together give the final 
difference equation for the differential equation (3.6,4). 
Considering equation (4.2.5), the difference equation for 
this can be written with the help of equation (4.2,2), as 


(I) = j ( 1 +1) ) + ^<K,( 2 - n ) 


rti. + ±3 ) 


in 


I m ■jjai.r.i. •§* I I« I 

^1 ^3 ^2 ^ 4 ^ 


hl+h3 I h^ hg 


\ hr 


2 + ^4 ) 


J 




(4.2.6) 


Where, 


T = JS J (M_ / 


"Y 


(4.2.7) 

Next, equation (4.2,5) is first re-written in 


the foim 




(II) = |oC (n - 1) 2 ) + 2 ^ 
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2 

Using the finite difference operator for ^ , this equa- 
tion reduces to 


(II) 


3 

2 


<<(n - 1) 



*%o. 



-f 
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hg" 



hg (hg4i^ (h^-^hgT T 4 


+ ^ 
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r 1 j. 1 1 * 1 

o ' hj_ hg hg h4 \ CY^'s^ 


-J^2_ %3 ^ .y4 ,r 

hg (hg+hp hg (hg+hp (h^Ht^y j-fo^ 


(4.2.9) 

Where, (at the node point i)- can be calculated from 
equation (4.2.7). It must be borne in mind that can 
never be negative. Combining equations (4.2.6), (4.2,7) 
and (4.2,9), the complete finite difference scheme for 
equation (3.6,4) comes out to be, 
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(1 +» ) + I »C (2 - n) j - 2 


ii- ' iXr-\ 

1 o 


11 + ±3] 

n r 




-C (n - 1) Y “'^1 ^2j^t_2Eo X + ^2 +-Xo' ' i 
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CO 
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7 ^>1 

+ ->2 
^2 ^ 



If, the mexh is assumed to be square, such that h^ = hg = 
^13 “ 5^4 “ h (say) then, equation (4,2.10) further reduces 
to 


(1 +D ) 


X (2 - n) 


- 4^ 
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Since, this finite difference equation is based on the 
values of ( or jC. ) at the nodal# points which have 

first to be calculated by finding out the derivatives of 
the stress function ^ , it is obvious that the solution 
can be obtained by an iterative process only. To start 
the solution a current set of values for will have to 
be generated. 

4,3 ELASTO -PLASTIC SOLLTIOH 

After developing the finite difference equations, 
both for the elastic and the plastic torsion, a procedure 
for obtaining the elasto-plastic ■ solution will be laid down. 
First we consider the criterion that has to be satisfied 
to obtain the required solution. Naturally, the elastic 
and plastic equations (3.6.3) and (3.6.4) should be satis- 
fied in their respective domains. A^;'. the boundary the 
stress function surface should be smooth so that the resul- 
tant shear stress is continuous at the elasto-plastic 
boundary. The computational procedure which has been used 
to obtain the solution can be summarised as follows: 

'(i) The region under consideration is covered 
by a rectangular grid. It is assumed that 
the surface passing over the nodal points 
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has nearly the same propert3.es as the actual 
stress surface, that is, its contours and 
gradients Kirill give the direction and the 
magnitude of the shear stresses respectively. 

(ii) The elastic torsion equation is first sol- 
ved over the entire section. To achieve 
this the difference equation is written 
corresponding to every node by using eqn, 
(4.2,2). In this way a set of algabraic 
simultaneous equations are obtained, which 
can be solved to give the values of ^ at 
the pivotal points. Since the number of 
simultaneous equations is generally more 
than a hundred, an iterative method has to 
be used for their solution. In this thesis, 
the Gauss-Seidel iteration method has been 
used, the details of which are given in 
appendix (a) . 

With the help of the se values of ^ at 
, X and hence the resul- 

E l-’l 

ss (T^ can be calculated 
byusing eqn. (4.2,7). 
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(iii) As a first approximation, we, can assume 
that the nodes at which the values of 
the non-dimensional stress {ov yi. ), 
exceeds one become plastic. Correspond- 
ing to such nodes the finite, difference 
equation is re-written using eqn. (4.2,10) 
and the previously calculated values of 

The difference equation correspond- 
ing to the other nodes remains unchanged 
since they still lie in the elastic 
region. 


(iv) This new set of simultaneous equations is 
solved for , After this? the values 
of and hence of ^ (or^ ) 

are calculated at every node again. 


(v) By using these values of the nodes 
are tested for being elastic or plastic 
by finding as to whether the values of 
are smaller or greater than one. 

Some new nodes may enter into the 
plastic region, whereas some of the nodes 
which were in the plastic region may 


become elastic 
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Therefore, the coefficient matrix of 
the simultaneous equations is again modified 
by re-writing the difference equation'.; for 
all ifhe nodes lying in the plastic zone by 
using equation (4,2.10). For the nodes which 
have become elastic (from plastic), the 
difference equation is re-written by apply- 
ing the difference equation (4.2,2). 

(vi) Steps (iv) and (v) are then repeated till 
the values of converge giving a stable 
stress distribution. In fact, the values 
of ^ and and hence of the elasto - 
plastic boundary will keep changing after 
every iteration. 

Depending on the number of iterations 
and of course on the mesh size any reason- 
able degree of accuracy can be achieved, 

(vii) Lastly, when the values of and have 
converged, the elasto-plastic boundary can 
be found by interpolation by using the 
criterion that the value of has to be 
one on it. Either a linear or a higher or- 
der interpolation formula can be used. 
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In going through the above steps (from (i) to 
(vii) , the elastic and the plastic equations (3.6.3) and 
(3.6,4) have been satisfied in their respective regions. 
With rqgard to the smoothness of the stress function <|> 
and the shear stress ^ , the method of finite differences 
itself implies that the derivative is continuous. No 
attempt has been made to justify this here. 



5 . BESULTS MD DlF;GUR,qiON 


5.1 appropriate MESH SIZE 

Since ti;ie numerica.! technique described in the 
previous section is based on finite difference approxima- 
tions, it becomes necessary -to decide an appropriate mesh 
size which will give a solution converging to the actual 
one . By carrying out the actual computation, it was found 
that the size of the suitable mesh depends on the shape 
of the cross-section. If in a direction, the gradient of 
the stress function (resultant shear stress) is expected 
to change rapidly, then a finer mesh must be used in that 
direction. This is necessary because the coefficients of 
the matrix in a row, corresponding to a node lying in the 

plastic region, are effected by the values of at its 

0 

neighbouring nodes raised to the (n-l) power. Therefoire, 
if the nodes are not sufficiently close, the diagonal 
nature of the matrix is not preserved and this leads to 
overflows on the computer. This type of difficulty was 
mainly experienced in obtaining the solution for a thick 
L-shaped section. 

For square sections the solution was obtained 
for different mesh sizes. Fig. 5 shows the variation of 
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the values of the stress function at the centre of the 
section with the step siz§. On the same graph the compu- 
tational time has also heen plotted for different mesh 
sizes. The stress function at the centre of the section 
has heen chosen as a representative value because it has 
the largest magnitude and is therefore most succeptible 
to change with the mesh size. 

From this graph it can he seen that after a cer- 
tain mesh size the solution converges to a value which 
does not change with a decrease in the mesh size . How- 
ever, the computational time increases very rapidly with 
the decrease in the step size. 

In our problem, therefore, the most appropriate 
mesh size for a square section, considering the accuracy 
as well as economy in time, is to take about 10 to 12 
meshes along the centre line OA, for the values of n less 

than about 20 . 

Similarly, for other sections also a proper mesh 
size may be chosen such that the method works and gives a 
converging solution with a minimum of computational time. 

5.2 CONVERGENCE TO ELASTIC IDEALLY PLASTIC SOLUTION; 

The method has been tested for convergence and 
accuracy, using the criterion that for large values of n 



54 


(or for large values of y\. for the piecewise linear appro- 
ximation), the solution should approach the elastic 
ideally plastic solution. That is, inside the plastic 
region, the resultant shear stresses should everywhere 

approach the yield stress = 1. 

^e 

The solution has been obtained for n = 30 with 

1.6. The maximum value of comes out to 
be 1.045 and 1.118 respectively. The contours of 
have been shown in Fig, (6) for 1.6. From this it 

can be seen that over most of the plastic region, the 
value of 5^ differs from 1 by less than about 8 percent. 

It is felt that if n is taken still larger, then the so- 
lution will approach the elastic- ideally plastic solution 
even more closely. However, we have not experimented with 
Values of n larger than 30 because of the large computa- 
tional time required (more than an hour on the I,B,M.7044), 

5.3 RESULTS FOR A SQUARE SECTION 

The proposed numerical technique has been first 
tried on a square section. This case has been investiga- 
ted in detail for various values of the angle of twist 
&/<c,^and the strain-hardening parameter n. This cross- 
section has also been solved for the piecewise linear 
stress-strain relation. 
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Figures (7) and (8) show the e.tasto-plastic 
boundary and the variation of, along the central line 
of the section respectively, for different values of • 

The effect of n on the elastoplastic boundary and on ^ 
has been shown in Figs. (9) and (lO). Fig. (ll) .shows the 
difference between the initial guess (obtained from the 
elastic solution) of the elasto-plastic boundary and the 
elasto-plastic boundary obtained from elasto-plastic so- 
lution. Some of the conclusions which can be drawn from 
these graphs are as follows: 

(1) The shapes of the elasto-plastic boundaries 
are sirailare for different values of the angle of twist 

5 3-^*^ they are almost parallel (see Fig. (7)). 

(2) Relative to its in3.tial position (corres- 
ponding to the first approximation based on the elastic 
solution), the elasto-plastic boundary does not shift 
along its entire length an one direction. It can be seen 
that over some of the portions near the centre line of 
the section, the elasto-plastic boundary shifts slightly 
into the elastic region inereasing the plastic area; 
whereas, over the remaining major portion it shifts into 
the plastic region towards the free boundary of the sec- 
tion thereby decreasing the plastic area- In Fig. (ll), 
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this shifting of 'the elasto-plastic boundary has been 
shown for the three values of/^*/feo= 1.2, 1.5 and 2.0. 

Fig. (9) shows that for n = 15 this change in 
the position of the elasto-plastic boundary is more than 
for n = 5, for all the three values of fact 

n = 5 gives a boundary which is very close to the Initial 
guess, obtained from the elastic solution. As n increases, 
this departure also becomes more and more . The nature of 
the shifting of the boundary remains the same as described 
earlier. Fig. (9) also suggests that small changes (by 
one or two) in the values of n will not alter the elasto- 
plastic boundary appreciably. Therefore, 

i) ForeO-l practical purposes, for the values of 
. ^/So slight ly mo re than one ( 1 . 0 

to 1.25 or 1,3), it is reasonable to assume 
that the region of plastic deformation is 
more or less the same as is obtained on the 
basis of the elastic solution only. This 
will be more so for materials for which the 
•walues of n are low. 

ii) For low values of n, the elasto-plastic boun- 
dary Is not sensitive to changes in n.^ 
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(3) However, Fig. (lO) shows that the resultant 
shear stress distribution within the plastic region al- 
ters app re ciabljr -with changes in the values of n, even 
though it does not effect the stress distribution in the 
elastic region. Hence the stress concentration will be 
different for different values of n. For example, for 
^/e-o= 1.5, ^ maximum (at the middle of the sides) will 

he 1.415 and 1.241 for n = 5 and n = 15 respectively. 

As the value of n decreases, the shear stress 
distribution approaches the elastic solution more closely. 
For large values of n, the solution converges towards 
the elastic-ideally plastic solution. 

Fig. (l2) shows the contours of the shear stre- 
sses over the cross-section. It can be seen that the shape 
of the contours do not seem to change even after plastic 
deformation has taken place. 

The non linear nature of the variation of the 
torque with the angle of twist ®/evhas been shown in the 
Fig. (l3). For a fixed 6 / 0 ^ , the torque remains almost 
constant for different n. 

The effect of changing the values of the para- 
meter o<_ in the stress-strain relation (3.4.5) has been 
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studied by obtaining the solution for - 1.5 and n = 15 
for 0i015, 0,020 and 0.025, Fig, (14) shows the values 
along the central line of the square section for 
all the three values of ^ , It can be seen that decreasing 
values of give higher stresses in the plastic region. 
Though the changes in the stresses are very small for these 
three cases, the values of the shear stresses and the 
elasto-plastic boundary will be affected quite appreciably 
if the permanent yield strain ^is taken to be very much 
different from 0.02 (say 0.002, as suggested by many 
authors). 

For the piecexirise linear work hardening approxi- 
mation the solution has been obtained for the values of 
the work hardening parameter 6 and l3.5 (i.e, corres- 

pondtog to E/E' = 5 and 10). 

Fig. (15) shows the elasto-plastic boundary, for 
®/( 9 q= 1.5, for 6 and n = 15, along with the initial 
guess of the boundary obtained from the elastic solution 
(i.e the contour along which, in the elastic solution, 

^ = 1. Note that for this guess the stress ^ will be 
higher than = 1 at many points). It can be seen that 
the shrinking of the elasto-plastic boundary relative to 
its first approximation based on elastic solution, towards 
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the free surface of the cross-section (decrease in plas- 
tic area) , is more for the piecewise linear work hardening 
approximation than, for the values of n around 15 in the 
polynomial work hardening approximation. 

Fig. (l6) shows the shear stress distribution 
along the central line of the cross-section. The maximum 
value of Q~ comes out to be as 1.087 and 1.044 for 
^=6 and 13.5 respectively. = 1,044, being very 

close to one shows that over the entire plastic region 
the resultant shearing stresses will be very nearly equal 
to the yield shear stress. Hence the values of ^greater 
than 10 ■'lay be assumed to give a solution which approaches 
the elastic-ideally plastic solution closely. This infe- 
rence is also consistant with the amount of shifting of 
the elastoplastic boundary for large values of n (See 
Fig. (16) for A == 6 and compare with n = 15), 

5.4 RESULTS FOR OTHER S.SGT IONS: 

The method has also been used for the solution of 
rectangular and thick L - sections. For these sections the _ 
plastic zone develops around more than one point over the 
cross-section. 
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For rectangular sections, solutions have heen 
obtained for the side's ratio of 1,2, 1.5, 2.0 and 5.0, 

For these four cases, Figs. (17, a),; (iV.h), (18) and (19) 
show the shape of the elasto-plastic boundary for =1.5 
and .n =15, It can be seen that the plastic region, which 
develops around the centre of the shorter side of the sec- 
tion, decreases with the increase in the side’s ratio. 

For a fixed the shear stress concentration is not 

much affected by changes in the side’s ratio. For exam- 
ple, the maximum value of for side’s ratios of 1,2 
and 5 is 1.209 and 1.193 respectively. Fig. (19) shows 
that for large values of the side’s ratio (b/a = 5) the 
elasto-plastic boundary, along the longer side of the 
section, remains almost parallel to the free surface, 

over the ma .1 or portion of its length. 

The solution has also been obtained for two types 
of thick L - sections. When the solution was first attemp- 
ted for the type of section shown in Fig. (20), by taking 
10 meshes along the vertical edge of the section, a solu- 
tion could not be obtained because of the ill condition- 
ing of the matrix after a few iterations. The procedure 
was then repeated after decreasing the me^ sisse till It 
gave a converging solution (for this case 14 meshes were 
taken along the vertical edge of the section); however 



61 


the solution was not checked for convergence with respect 
to the mesh size. While obtaining the solutions for other 
sections, it was found that if the solution converges for 
a particular grid spacing then it is close to the actual 
solution 5 othervn.se the solution does not converge at all; 
in fact for a square section the elastic-plastic boundary 
is hardly affected. The other type of L - section, shown 
in Fig. (21), has been solved for two values of n (n = 5 
and n = 10). For this section, when 15 meshes were taken 
along the vertical edge of the section the matrix was 
much smaller than for the previous L ~ section and less 
computer time was required. 

Therefore, it may be concluded that, for diffe- 
rent shapes of L - sections, different types of grid 
spacings are required to obtain a converging solution. 

As a suggestion for further work, it maybe of use to 
take a finer mesh around the point of high stress concen- 
tration where the stresses vary comparatively more abrupt- 
ly; over the remaining portions of the section where the 
stress variation is slow, a comparatively coarser mesh 
will not affect the solution. 
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Figs. ( 20 ) and ( 21 ) show the shape of the plas- 
tic rogions developed over the cross-sections of thick 
L - shaped sections, obtained by cutting out a square 
from a square, for 0 / 0 ^ = 3 , 

As expected, the elasto-plastic boundary near the 
point of maximum stress concentration suffers more change 
in comparision to the other regions of the cross-section 
undergoing plastic deformation. The shear stress concen- 
tration is affected very much by the elasto-plastic de- 
formation. Due to this the stress distribution in the 
plastic region, around the point of maximum stress concen- 
tration, becomes more iinlform than in the purely elastic 
case. The maximum value of Tf- for the L - section shown 
in Fig. ( 21 ) comes out to be 1.649 only for '^/©o “ ^ 
and n = 10 . 

5.5 CONCLUDING REMARKS 

The numerical method which has been described in 
section 4, for obtaining the solution of elasto-plastic 
torsion problems, can be used for any other cross-section. 
The success of the method in obtaining the solution for 
the L - section shows that it is quite capable of giving 
the solution for sections having very high stress concen- 
trations at some points on its boundary. The only care 
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that needs to be taken is to decrease the mesh size near 
such points of stress concentration. The greatest ad- 
vantage of this technique is that it automatically takes 
care of all the plastic zones developed over the section 
and gives a solution satisfying the elastic and plastic 
torsion equations in their respective regions. 

The method can also be used for other work - 
hardening laws in which is some other polynomial func- 
tion of ^ i.e. -€ = f ( It needs to be further 

tried for multiply connected and hollow sections. 
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Diagonal systems of simultaneous linear alge- 
braic equations have the fundamental property of being 
efficiently solvable by iterative methods. Among these, 
the Gauss - Seidel iterative method has the advantage 
of being simple to use. To apply this method, each equa- 
tion of the system is first solved for the unknown with 
largest coefficient; so that the equation may be written 
as, 


^1 ^12 ^2 ^13 ^3 


a,n X + k 
i n 1 


X, 


^23 ^3 


. agn x„ + kg 


* # • « 


• « « « * 




n 


\ l 1 ^1 ^ n 2 ^2 


^n,n-l ^n-1 

( 0 ) 


For using the iterative process, initial values Xj 
for the unknowns are first assumed for starting the ite- 
rative process. New values of x^^^^ are then obtained 

from the left hand side of Eqn. (A.l). These new values 
( 1 ) 


X.. ' are substituted in the right hand side of the 
3 ( 2 ) 

equations to give the second approximation x^ . This 

process is continued untill is equal to x^^^ 

to within the required accuracy. The are then 

the roots of the system. 
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In Gauss's method each approximation of a 
root is obtained hy means of a single m.achine operation^ 
and errors do not impair the convergence of the process 
since they are equivalent to a new set of starting 
values. The method, if convergent, converges no matter 
what starting values are chosen. In this method, at 
any stage, the best available approximation of the un- 
known is used in the right hand side of eqn. (A.l). 
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APPENDIX B 

gRIg£.igJCRIPTI ON OF T HE FORTRA N PRO GR At4 FO B. T HE SOLUTION 

OF BLACTO- PL ACTIO 'PROBLEM 


The niimerical procediire outlined previously 
for solving the elasto-plastic torsion problem was pro- 
gra«d for solution on the I.B.M. 7044 FORTRM IV. 


The correspondence between the F0RTRi\N sym- 
bols Used in the program and the mathematical symbols 
employed previously is shown in the following table. 


FORTRM 

MATHMATICAL 

f FORTRAN 

mathematical 

SYMBOL 

SYMBOL OR 

SYMBOL 

SYMBOL OR 


DESCRIPTION 


DESCRIPTION 

A 

Matrix coefficients 

THETA 

®/eo 

o 




X 


YRE 

(200) h 

DX 

DY 


YRP 

YIELD 

(E © ) h^ 

(yield stress 

Y i 

. : ' ", 1 

i 

Constant term b^ 
in the R*H.S. ■ 
of the difference 
equations 

H 

N,M 

in tension) 

h (mesh size) 

No of finite 

TOIaF 

q (resultant 


meshes taken 


shear stress) 


along or 
axis. ^ 

CHI i 

0 % 

NP 

" ' 1 

No. of columns 
of the banded 


' 


matrix. 
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FOKERAIT " ' 
SMOL 

1 mathem^ical i 

SYMBOL OR j 

DESCRIPTION i 

FORTRAN 

S3!MB0L 

! 

! 

j mathematical 

1 SYMBOL OR 

1 DESCRIPTION 

ITUM 

j 

! 

No. of iteration 

1 

i 

XOLD 

1 

i 

1 

1 Values of 
in the previ- 
ous iteration 

ITURNL 

Maximum or last 

1 1 

AK 

No . of iteration 

1 ! 
CHIOLD 

Values of cV 

( 1 +1) ) 

Or 

CHIl 

in the pre- 
vious itera- 

AIFA : 

X 


tion. 

YAN 

n 

KNl 

Number of lo- 
west nodal 

TORQUE 

T 

KN2 

i 

point along a 
vertical 
grid line. 

Number of top 
nodal point 
along a ver- 
tical grid 

1 line 

I 



< ■ i 



The following remarks are intended to help with 
the Tinder standing of the program. 

(1) Before starring the program, a net of grid lines is 
spread over the cross- section and the nodal points 
are numbered in the manner shown in Fig. (22). 

(2) Subroutine ELMAT is used to write the coefficient 
matrix of the difference equations corresponding to 
the elastic torsion equation (3.6.3). 



While solving such partial differential euqa- 
tions by finite difference method, the coefficient 
matrix is bah.ded. Hence, the coefficients of the 
matrix have been stored in the memory by re- stack- 
ing the band in a vertical manner. This is shown 
in the figure below. This procedure reduces the 
memory re qui rement s . 
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(3) Subroutine GAUSS has been written to solve the above 
set of simultaneous equations. 

(4) Subroutine STRESS is used to calculate the effective 
stress ^ at all the nodal points lying inside the 
region and at the boundary of the cross-section. 

For calculating sX the boundary points, where- 
ever necessary, partial derivatives of ^ have been 
computed by using one sided three or four point 
formulae . 

(5) In going from one iteration to the next and for cal- 
culating the stresses, instead of talcing the new 
values of ^ , the old values of $ (from the pre- 
vious iteration) are changed by a fraction of the 
difference between the current and the previous 
values. This is the reason for storing the values 
of ^ of the previous iteration in the memory. 

(6) Subroutine PLAST has been used to modify the coeffi- 
cient matrix. Every node is tested, one by, one, to 
check as to whether it lies in the ‘.•elastic or 3n the 
plastic region. ‘Whenever, a node point shifts from 
the elastic to the plastic region or from the plas- 
tic to elastic region^ the corresponding difference 
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equation (coefficient of the matrix in that row) is 
re-written by using the difference equation (4.2.11) 
or (4,2,3) depending on whether it shifts from the 
elastic to the plastic region or vice-versa. 

(7) Subroutine TORQ has been used for calculating the 
torque (T). For this the trapezoidal rule of double 
integration has been used. 

( 8 ) The last subroutine IRTERP carries out the interpo- 

lation of the elasto-plastic boundary along every 
vertical as well as horizontal grid line, using the 
criterion that 7 ^ = 1 at the elasto-plastic 

boundary. 

FLOW DIAGRAil 

A flow diagram of the procedure employed is 

shown on the next page. 
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FLOW diagram 



Start 


Read 


%Oy ITURC^L, IVIJl, IVIJ2, IVIJ3 | 


call elmat 

(Generating the finite difference equa~ 
tions at every node for elastic torsion.)! 


call gauss 

(Solution of the difference equa- 
tion using Gauss-Seidel method.) ' 



ITURN>0 



Fract 

( ^old-f ) 


call STRESS 

(Calculation of the_re suiting shear- 
ing stress q and at every node.) 



gall PLASr 

(Testing every node whether it lies in the elastic 
or plastic region and correspondingly modifying 
the difference equations nsing equation(4.2.1l) 

(difference equation for plastic torsion)), 

' ’ ' ' ' ' ' ’■ ^ 

V 


Repeat the step 3, 4 & 5 | 

till the values of con- 
ve rge . ^ 

t 

alculate the resultant shearing stresses. 


X 


Print # , err 

1 1 

'f 

1 — 

gall TORQ 

(Calculation of Torque 


gall IHTERP j 

(Interpolation of the elasto-plastic boundary i 
along every horizontal and vertical grid line.) I 
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ELAST-PLASTIC SOLUTION OF A SQUARE SECTION 

ITURN=Q 

ITIN=C 

READSllt TliETA, ITURNL, IVIJl, IVIJ2,IVIJ3 
PRINT BlltTHETA. ITURNL» IVIJl, IVI J2, IVI J3 
311 FCRMAT{5X,FA. 1,4( IX, 12 ) ) / 

N=1C 

V:=N 

H=l./ (24»V1 

NL-={N«*2-N)/2+N 

IP=N-1 

NP-2*IP+1 

IVIJP=IVIJ1+1 

YRE=1C0. 

OIRENSION A( 55, 19 ) , X( 55 ) , Y( 55) , DX{65) ♦0Y(65) ,CHI (65) 

CIRENSION TCW(65),XOLD(55 ),CHIl{65) rVERTI 11) ,HORZ(ll) 

CO 1 1=1, NL 
CC 1 J=1,NP 

1 At I,J)=0. ’ 

■DG" 2 1 = 1, NL ■ 

2 x<i) = iooo.- 

DC- 3 1 = 1, NL 

3 n I ) = YRE 

CALL ELMAT(N,NL,AtNP,IP) 

351 CALL GAUSSIX, Y, A,NL, IP,NP,N, ITURN) 

IFi( ITLRN.EC. I.QR.ITLRN.GT. DGQ TO III 
353' PRINT 342 

342 FCRMAT{//5CX, 13HVALUES OF PHI) 

PRINT 343, tX( I ) , I = ltNL ) 

343 FCRMAT(/X,8EL6.6 ) 

CALL STRESS ( X, nX , DY ,NL , H, N , TOW,CHI , YRE ) 

NL1=NL+I 

NL2=NL+N 

IF t IT INi EQ.O) GO TO 600 
GO TO 660 

600 ■ VRE=YRE/CHI { NL 1 ) *THETA 
YRP=5l .25»YRE 
DC 61C I = 1,NL ' 

610 Y(I)=VRE 

DC 605 I=1,NL 

605 XU ) = X ( I) /CH n NL I ) *THETA 
IT IN= IT'I N+ 1 , 

1^/3 ' ' ' T'"'' ' 3''2 S '''''''' '''"' ' ' ■' ' ' ■ ' Vvv^;, 

325 F0RMATt//4QX,49HVALUES OF RESULTANT SHEAR STfl^^ES AT THE BOUNDRY) 
PRINT32T,ITDWt 1) , I=NL1,NL2) 

327 FQRMAT(/1X,16F8.1) 

PRINT 323 
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ELAST~PLASTIC SOLUTION OF A SQUARE SECTION 


311 


1 

2 

3 

351 

353 

342 

343 


600 

610 

605 

660 

325 

33i 


« If «#-*!## If ^ It "If If If II ###### ##||« I 

ITURN=Q 

rTIN=C 

REA0311, TIlETAt ITURNL , IVIJl, IVIJ2,IVIJ3 
PRINT 311tTHETA, ITURNL, IVIJl, IVI J2, IVIJ3 
FCRMAT<5X,F4.1t4{lX, I2) ) / 

N=10 
V = N 

F=l./ (2;»Vl 
Kl;={N«*2“N)/2+N 
IP=N-1 • 

NP=2*IP+1 : 

iyijp=iviji+i 

YRE=1C0; 

OIRENSION AlSSr IR ) » X155) rY(55) ,DX(65) ,0Y(65> ,CHn65) 

DIMENSION T0W{65),X0LD(55),CHIl(65},VERT(ll)THaRZ{ll) 

CO ,1 I=1,NL 
DC 1 J=1,NP 

A( It J )=0. ' 

■DC 2 1 = 1, NL 
X-( I ) = 1000. 

DC- 3 1 = 1, NL 
Y1 I)=YRE 

CALL ELMATINtNL, A,NP, IP) 

CALL GAUSStX,Y,A,NL, IP,NP,N,ITURN) 

IFi( ITLRN.EG. l.OR. ITLRN.GT. DGO TO 111 
PRINT 342 

FCRMAT(//5CX, 13HVALUES OF PHI) 

PRINT 343, IX( I ), I = 1,NL) 

FGRMAT{/X,8E16.6) 

CALL STRESS(X,nx,DY,NL,H,N,TOW,CHI,YRe) 

KLl=NL+l 

Ni2=NL+N 

TF{ ITINiEQ.O) GO TO 6C0 
GO TO 660 

YRE = YRE/CHnNLI )*THETA 
YRP-1.25»YRE 
DC 61C 1=1, NL 
YI I )=YRE 

DO 605 1 = 1, NL , ' 

X.n) = X(n/CHnNLU»THETA 

IiTIN=ITiN+l 

■GO TO 35 3 

FCRMAT(//4QX,4RHVALUES OF RESULTANT SHEAR STRESSES AT THE 80UNDRY) 
PRINT327,(T0W( n,r=NLl,NL2) 

F0RMAT(/1X,16F8.1) 

PRINT 323 



XI 


323 FGRMAT{//45X,34HVALUES OF RESULTANT SHEAR STRESSES) 
PRINT327, (TOW{ I ), I=1,NL) 

PR INT 33) 

PRINT335,(CHI( I),I=NL1,NL2) 

PRINT 333 

PRINT335i(GHI{I),I=l,NL) 

331 F0RMAT{//4GX,27HVALUES OF SIGMA AT BOUNDARY) 

333 FGRMAT(//50Xt ISHVALUES OF SIGMA) 

335 FCRMAT{2Xtl6F8.3) 

IF( ITURN.EG.0.CR.ITURN.GT.IVIJ2) GO TO 777 
GC TO 775 

777 CALL TORQtX, H, TORQUE, Nt IP ) 

CALL INTERP(CHI,VERT,HORZ,IP,N,NL) 

775 DC 35 C I=1,NL 
350 XCLD(I)=X(I) 

IF( ITURN.EQ.OIGO TO 361 
GC TC 367 

361 DC 363 I=l,NLv 
363 CUK I ) = CHI ( I ) 

367 ITURN=ITURN+1 

IFU-TURN-ITURNL) 354, 354,357 
354 PRINT 358 

358 FQRMAT(1X,128{ IH*) ) 

PRINT 359, ITURN 

359 FCRMAT(//6QX,6HITURN=,I3) 

CALL PLASTI IP,N,NL,NP,A,CHI,Y,CHIi,YRE,YRP) 

GO TO 351 

111 fF{ ITURN. GT.IVIJU GO TO 113 
IF< ITURN. GT. IVIJ2)G0 TO 115‘ 

121 IFt ITURN. GT.IVIJ3) GO TO 119 
FRACT=0.2 
GC TO 117 
119 FRA€T-0i5 
GO TO 117 
115 FRA:CT=Q.8 . 

GO' TO 117' ■ 

113 FRACT=l; 

11? CO 352 I=1*NL. 

352 X,( I)=XOLD(I) + FRACT*(X(I)-XOLD{I) ) 

GO TO :.35 3 ^ 

357 CALL TORQIX, H, TORQUE, N, IP) 

CALL INTERR(CHI,VERT,HORZ,IP,N,NL) 

■ : STOP 

ilBFTC SU81 

SUBROUTINE BLMAT ( N ,NL , A , NP , I P) 

DIMENSION A(NL,NP) 

A{1,1)=4. 

Ai(l,2)=-4. 

DC 1 1=2, IP 
A4 I, I-l) =-X. 

AI I, I )=4. 

A:{ I, I+l)=~l. 
n=IP + I 
Ai{ I, 1 1)^»-2. 



XII 


A(N,N )= 4 . 

MI = N-HP 
A{N,NII)=~ 2 . 

KE = 0 
KNl = iSl + l 
IRRO=C • 

L=N • 

NS=NL-IP 

Nl=N -2 

CC 9 I= 1 ,N 1 

IR=I +1 

TRR=I 

IRRN= IRRG+IRR 
KN 2 = N'«IR-IRRN 
■ ^ 1 = 1-1 

CG 2 J=KNltKN2 
N?f'L = N-L 
NPL=N4L 

If { J.GT.NSJGO TO 4 
A"{:j,N)= 4 . 

I..F( J.EQ.KNIIGQ TO 3 
ATJtNf'Ll =-l. 

A{J,NFL- 1 )=- 1 . 

A«J,N -13 =-l. 
inj.EC.KN 2 >GO TC 2 
A ( J,N + 13 =-l. 

GC TO 2 

3 AIJrNRLT=-2. ' 

A ( J,N + 13 =- 2 . 

■G'G TC 2 ■ 

' 4 ^ • 

NRNE=h+NE 
NMLP=N-L+NE 
RPLP=N+L+NE 
A( JtNPNE)= 4 . 

IFiJ.EC.KNl) GC TO 5 
A(J,N^'LP3=-l. 

A{ J»NFLP- 11 =- 1 . 

A(J,NPNE- 11 -~ 1 . 

TP^J.EQ4KN2)Ga TO 2 
A (JtNPNE + ll=-l- 
: ' GO TO ' 2 \ ■ 

5 AtUfNPLP )=-2. 

A(.J,NPN6 + U--2. 

2 CONTINUE 

T ::v. : - 

v::::9::"::::/^V;gO^NTTRUE-" 

SUfiROLlTlNE GAUSS(>!»Y,A»NLfIP»NPfNf ITURN) 
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XIII 


CINEKSION X{NL),Y(fs]L),A(NL»NP) 

IT = 1 
ITL=2CQ 

. tf(ITURN.GT.I) go to 21 
GO TO 20 
21 rTL=75 
20 DC 60 1=1, N 
P = Y( I ) 

DC 50 J=1,NP 
If ( I~J}40,5C,4C 
40 P=P-A{I, J)*X( Jl . 

50 CONTINUE 

X( n = F/A( I, I) 

60 CGNT INUE 
INC=1 
MT=N+1 
NTR=NL-IP 
■ DC 70 I=NT,NTR 
P = Y ( I ) 

DC 65 J= 1,NR 
INi\=INC+J 

IF( I-{INC+J) )68,65,68 
68 R = P-Aa, J)»X(IKN} 

65 CONTINUE 

X ( I ) = P/A( UIP+1) 

INC=INC+1 
70 CONTINUE 
INC=INC-1 
NNP = NL-IP-F1 
NH=N+1 

DC 90 I=NNP,NL 
P = Y( I ) 

DC 80 J=1,NP 
NPP=NL-NP+J 
IF(I-(INC+J))85,80,85 
85 P=P-A{I, J)*X{NPP) 

80 CONTINUE 

X( n = P/A( I,NH) 

NH=NH+1 
90 CONTINUE 
I T = I T 1 

30 . RETURN 

• - -END • ^ 

£IBFTC SUB3 

SUBROUTINE STRESS! X,DX,DY,NL, H,N, TOW, CHI ,YRE) 
D I#'ENSION Xt55) ,OX(65) ,0Y{65) ,T0W165) , CHI (65) 

CALCULATION OF DERIVATIVES IN Y-DIRECTION 



u u u 


XI? 


: m=i 

JRRN=JRRO+JRR 

LNEW2=N*JR-JRR'N 

DC 2 J=LNEW1,LNEW2 

IF( J.EQ; LNEWl IGO TO 3 

IF{ J.EQ. LNEW2)G0 TO 4 

DY{'J) = (X{ J + 1) )/{2.*H) 

GO TO 2 

3 IF(J.EQ.NLMT)GO TO 5 
IFU.EQ.DGC TO 6 

DY(J)={3.-»X{J)-4.»X( J+1)+XT J + 2) )/(2.*H} 

GO TO 2 

4 CY(J)=(X(J-1)-C.)/(2.*H) 

GO TO 2 

5 CY( J)={3.*X(J}-4.*X(J+1) )/{2.*H) 

GO TO 2 

6 CY{1)=0^ 

2 CONTINUE 

LNEW1=LNEW2+1 
JRRO=JRRN 
1 CONTINUE 

DYtNL ) = ( 3.*X(NL)-4.*X(NL-l)+X{NL-3) )/ {2.*H) 

NEX1=NL+1 

NEX2=NL+N-I 

K=N 

J = 0 

DC 15 I=NEX1,NEX2 
J=J + K 

DYI l)=-(4.*X{J)-X( J-1) )/{2.*H) 

15 K=K-1 

NLRN=NL+N 

OY.<NLFN)=~(YRE + X{NL) )/{2.*H) 

CALCULATION OF DERIVATIVES IN X-DIRECTION 

23 CX(1)=Q. 

DC 25 I=NEX1,NLPN 
25 CX:U)=a. 

LNEW1=N+1 

' Jft«0=C ^ :■ ■ : 

DC 27 I=ltNMT 

v-:-: J:RR=:c : 

JR#N=JRR:0+:J^RR - 
LNEW2=N* JR-JRRN 

IF( J.EClLNEk^DGO TO 31 

OXt J) = (X( JN!LS)-X{ JPLS-1 ) )/t2.*HJ 



XT 


31 DX(:J)=DY{ J1 
29 CCNTINUE 

LKEW1=LNEW2+1 
J{?RO = JRRN 
27 CCNTINUE 

DX(NL)=DY{NL5 

C CALCULATION OF SIGNA AT DIFFERENT POINTS 

YIELD=35CCG. 

DC 33 1=1,NLPN 

33 CFK 15=3. f*G.5»{DXU)«*2 + DY{I)**2)**0, 5/YIELD 
CC 35 I=1,NLPN 

35 TCW(I5=(DXtI)**2+DY{I5**2)**0.5/l. 

RETURN 

END 

£18 FTC SUB4 

SUBROUTINE PLAST UP, N,NL,NP,A,CH I, Y, CHI 1»YRE» YRP 5 
CIKENSION A{NL,NP),CHI{65),Y(NL),CHI1(65) 

AK=1.25 

YAN=5. 

ALFA=C.02 
CO 1 1 = 1, IP 
ri=rp+i 

IF(CHI(n-l.)ll,ll,52 
11 IFICHIK I ) .LT.l. )GO TO 1 

Y ( n =YRE 

A{ I, I-l)=-l. 

AU,I+1)=-1. 

AU,II5=-2. 

A.{ I, I 5=4. 

^GC TO 1 

52 CONSl=AK+1.5*ALFA*( 2.-YAN)*(CHI (I )**( YAN-1. ) ) 
CONS2=1.5*ALFA*(YAN-l. ) * ( CH I { I ) ** { YAN-2 . ) 5 

Y l I ) = YRP 

A { I, I~1)=-U.*C0NS1+C0NS2*(CHI ( I 5+CHI I I-l) 5/2. ) 

A( I, I-H) =-{;l.*C0NSl+C0NS2*(CHI(I)+CHI (I + l) )/2. ) 

AU, I !)=-( l.*CCNSl+C0NS2*(CHI(n+CHI( II) 5/2. )*2. 

A{ 1,1 )=-( At I, I-l)+A( I, I + 1)+A(I,II ) ) 

1 CCNTINUE 

CGNSl=AK + 1.5^fALFA*(2--YAN)*{CHI (N)«*(YAN-1. 5 ) 
CCNS2=i; 5*ALFA* ( YAN- 1.5* (CHI {N)**{ YAN-2. ) 5 
Ai(N,N-l) =-t l.*C0NSl+C0NS2*tCHI (NJ+CHI (N-l) 5/2. 5 
NPIP = N+'IP 

A<N,NFIP)=-I l.«C0NSl+C0NS2*(CHI(N)+CHI{NPIP) 5/2. 5*2. 
EXTRA =-( l.*CONSl+CONS2*{CHI (NI+CHUNL+I) 3/2. 5 
A(N,N )=-{ A{N,N-1>+A(N,NPIP5+EXTRA) 

Y;(N)=YRP 

IRRO=C 

/U-y- : ;^-:y:>:::;yy 

,:;:'^y:y:c:^^yy:DtikMy:iMi;YN:i:'y-::y^^^::‘y;:'^ 



: ^ : vV Xn 

L-L-1 

NN'L=N-L 

KPL=N+L 

KLPIR=NL + IR 

CC 2 J=KN1,KN2 

JN1=J-L 

JPL=J+L 

rF(CHI(J )-l. ) 12,53,53 

12 ' IF{CH IK J ) .LT. 1. )G0 TO 2 

Y( J)=YRE 

IFtJ.GT.NS 3.G0 TO 14 
)=4. 

IFTJ.EQiKNl) GO TO 13 
J,N-1} =-l. 

4tJ,NPL)=-l. 

AU,NPL-13=~l. 

IFt J.EG^KN2)GC TO 2 
/i{4,N + 13 =-l. 

GO TO 2 

13 A{J,N+l)=-2. 

A(U,NPL)=-2. 

PC 102- 

14 NE=J-NS 
NPNE=N+NE 
NW1P=N-L+NE 
NFLP=N-KL+NE 
AilJ,NPNE)=4. 

IF^ J.EGlKNll GO TO 15 
At J,NFNe-l)=-l. 

A.(J,NFLP )=-l. 

A{J,NFLP-11=-1- 
IF(J.EG.KN2)G0 TO 2 
A,(J,NPNE + 11--1. 

GC TO 2 

15 A{J,NFNE + l).=-2- 
A.(J,NFLP )=-2. 

60 TO 2 

53 CONS i = AK + 1 .5*ALFA* ( 2 .- VAN ) * tCHl t J ) ■»* t YAN- 1 . ) I 
CCNS2 * U 5 * ALT- A* ( YAN- 1 . ) * ( CHI ( J1 ** t YAN-2. ) 1 
Y1J}=YRP ; ^ / 

IFl J.GT.NSIGO TO 4 
IF^tU.EQiKNl} GO TO 19 

A(.J,N-l)=“tl.*C0NSl + C0NS2*{CHI{ J)+CHI { J~l) )/2. ) 
AtJ,NHL)=-ti.*C0NSl+C0NS2*tCHKJ)+CHI( JML) )/ 2 . ) 

A t J,NPL-l)=-( l.*C0NSl + C0NS2*tCHH J)+CHI t JPL~1) )/2.) 
IFjtJ.EG-KN2}, GC TO 3 

At J,N + l)=-til.*C0NSl + C0NS2*tCHIt J)+CHI t J + 1) 3/2. ) 

A t J,N)=-(AIJ,N-1 )+At J,N+1)+At J,NML)+At J,NPL-1) ) 

3 EXTRA=-t l.«C0NSl+C0NS2*tCHIt J)+CHItNLPIR) )/2.) 

AtJ,N)=-t AtJ,N-l)+At J,NML)+At J,NPL-1}+EXTRA) 

19 AlJ,N + U=-t.l.*C0NSl+C0NS2«tCHIt JKCHIt J+l) )/2. )*2. 

At.J,NML3 =-Cl«*C0NSl + C0NS2*tCHI{ J)+CHI t JML) ) /2. )*2. 



mil 


GC TO 2. 

4 NE'-J-KS 

NPKE = N+NE : : ' . 

NMLP=N-L+Ne 

NPLP=K+L+Ne 

■ IF(J.EQ.KNl) GC TO 21 

A,( JtNFNE- !).=-{ l.»CONSl + CONS2*{ CHI U )+CH I (J-l) )/2. ) 

4( J,Nf'LP3=-( l.»CONSl+CONS2»(CHU J)+CHn JML) )/2. 5 
A{J,NPLP-ll=-{ l.*C0NSl+C0NS2*{CHI { J)+CHI { JPL-1 ) )/2-) 
IF{J.EQ.KN2) GC TO 5 

^ U,NFNE + l)=-( l.*CONSl + CONS2»(CHn J3+CHI ( J+1) )/2. } 
ii{ J,NPNE)=-{ A( J,NPNE-I)+A( J,NPNE+13 +A{ J,NPLP-1)+A{ JtNMLPT) 
: GO TO 2 ’ 

5 EXTRA=-{ l.*CON'Sl+CONS2*(CHI { J >+CHI (NLPIR) }/2- ) 

AT JtNPNE)=-£ AC J,NPNE-1)+A{ JtNPLP-1 3+A( J,NMLPJ+EXTRA3 

GO ro 2 

21 AtJtNPLP)--Cl.*CONSl+CONS2*CCHl{JJ+CHI{JML5)/2.)*2. 

AC J,NFNE + 1).=~{ l.*CONSl+CONS2^^(CHl{ J)+CHI( J + 1) )/2. )*2. 

AC J,NFNe)=-{ AC J,NMLP )+A{ J,NPNE+U ) 

2 CONTINUE 

KM = KN2+1 
IPRO=IRRN 
•9, CONTINUE 

RETURN 
ENC ■ : 

£IBFTC SUBS 

SUBROUTINE TORQCX,HtTORQUE,N, IP) 

DIMENSION XC55) 

TGR0UE=0.5*XC 1)+ 

DC 1 1=3, N : 

1 TCR«UE=T0RQUe+2.*XCl) 

KN1=N+1 ■ 

IRRD = C 
DC 2 1=1, IP 

IR=I+'l : 

IRRN=IRRC+I . 

KN2=N*IR-IRRN 
" DC 3 J=KN1,KN2 

IFCU.EQ.KNllGO TO 4 

: : knR1=kni+i : : - : 

: : iFiu:.EiiK:NPT)coy.TO-:-:B : : 

: T0BfiUE»TGRSUET^'4.*^^^^^^^ .. ,9. . : > 

GO ' TO’ :3: '• /V v 

4 TGRQUE=TCRCUE+2.*XC J) 

S:- ;-vTCR«UB=TORO:UE+4.;*X CU 

■>:;3;;’9-’’:cont : 

’ • ::’’9 : / : 

;:-:2::’'’’’9rcoNisiNUE:: 9 ', 

TORQUE=TORQUE*H**2/4.*8. 

5'’^’-::;::’;’:'"^;®R1’NT’®3B:0WC:R0UE’^^ 

370 F0RMATC///45X,*TDRQUE=*,F12.3) 
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£IBFTC SUB6 

SIjBRQUTINE INTERP(CHI, VERT.HORZt IP,N,NL) 
DI^'ENSinN CHI (65) tVERT C 11) ,HORZ{ 11) 

TMERFOLATIGN GN VERTICAL LINES 

KM = 1 
JRR0=C 
DC 1 1=1 , IP 
jRR=I-l 
JRRN= JRRC+JRR 
KK2=N«I-JRRN 
NLPr=KL+I 

IFICHKNLPI )«LT. 1, )G0 TO 2 
IF(CHI(NLPI).EQ. l.)GQ TO 3 
J=KN2 
V = Q. 

4 TFtCHUJI.LT.l.IGC TG 5 

IF(CHI(J).EQ. l.)GO TC 6 

J=J-1 

V=V+1. 

GC TC 4 

6 VERT( I)=0. 

^GC TO 2 ■ : 

3 VERT( n = 0. 

GO TO 2 

5 IF{J-EQiKN2).GO TO 21 

^ERT( n=V+tCHI{ J+l)-l. )/(CHI(J+l)-CHI( J) ) 

''gl] ■■ ff Q"' ,"2' ■ 

21 .VERT{ I) = V+CCHI{NLPn-l. )/(CHI(NLPn-CHI( J) ) 
2 KN1=KN2+1 

JRRC=JRRN 
1 CONTINUE 

INTERPOLATION ON HORIZONTAL LINES 

NP=N+1 
OO 7 1*2, NR 
IFtl-EQ.KPlGO TO 13 

’i- j ' ^ ■ ■■ . y y ^ 

IF.{CHI(I).LT.l.O.OR.CHI{I).EQ-l-O)0O TO 7 

r-v j*a.+L ■: . 

IFICHKJ )-LT.l. )GO TO 9 

ifcchhj).eg.i-)go to b 

i:.- V ^ 

9 HORZ( n=H+tCHl{J8)-l-)/(CHi(4B)-CHI{J) > 

8 HCRZt I) = H+l. 



XU 


Hi=0. 

15 IF(CHI{J).LT.l. )G0 TO 17 

IF(CHI{ J ).EQ. 1. )G0 TO 19 ' 

J =U + 1 
H=H+1. 

:GG TO \15 

19 HCRZ( I) = H41. 

GO TO 7 

17 HCRZ( n = H+lCHI ( J-D-l. ) /{CHIl J-1)-CHI ( JH 

7 CONTINUE 
PRINT^25 

425 F0RMAT{//35X,*VERTICAL DISTANCES OF THE BOUNDRY*) 
PRINTASn, { VERT( 1) , 1 = 1, IP ) 

430 FCRMATC/ 1X,15F8.3J 
PRINT435 

435 FCRMAT(//36Xt*HORIZONTAL DISTANCES OF THE BOUNDRY#) 
PRINT439, (HORZt I ) t I=2,NP) 

439 FCRMAT(/ lX»16Fa.3) 

RETURN 

END 


SENTRY 
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(Fig. 4b) PfECE-WISE LiNEAR STRESS -STRAIN CURVE. 
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6) CONTOURS OF SHEAR STRESS Oi IN THE 
PLASTIC REGION. 












(Fig. 9) ELA5T0-PLAST1C BOUNDARIES OF A SQUARE 
SECTION FOR DIFFERENT VALUES OF n, 



X/OA 


(Fig. 10) NON-DIMENSIONAL STRESS, Ol ALONG THE 
CENTRE LINE OA FOR DIFFERENT VALUES 
OF n AND O/Go =1-5. 






--- ELASTIC SOLUTION 
ELASTO-PLASTIC SOLUTION 


(Fig. II.) ELASTO-PLASTIC BOUNDARY OF A SQUARE SECTION 






TORQUE 



1-75 


TWIST' 


1-0 

ANGLE OF 


(Fig. 13) % 




(Fig. 14) m ALONG THE CENTRAL LINE OA FOR DIFFERENT 
: : VALUES OF dC.;.' : 
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o- POLYNOMIAL WORKHARDENING (n = 15) 


ELASTIC SOLUTION (1st. APPROXIMATION) 

(Fig. 15) ELASTO-PLASTIC BOUNDARY FOR PIECEWISE 
LINEAR STRESS STRAIN CURVE. 







i . iELASTO-PLASTIC BOUNDARY FOR RECTANGULAR 
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(Fig, 19) ELASTO-PLASTIC BOUNDARIES IN A RECTANGULAR 
SECTION b/a = 5. 






(Fig. 21) ELASTO-PLASTIC BOUNDARIES IN A L-SECTION CUT OUT OF 
A SQUARE. 













